Disorder-enhanced phase coherence in 
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The consequences of disorder on interacting bosons trapped in optical lattices are investigated by 
quantum Monte Carlo simulations. At small to moderate strengths of potential disorder a unique 
effect is observed: if there is a Mott plateau at the center of the trap in the clean limit, phase 
coherence increases as a result of disorder. The localization effects due to correlation and disorder 
compete against each other, resulting in a partial derealization of the particles in the Mott region, 
which in turn leads to increased phase coherence. In the absence of a Mott plateau, this effect is 
absent. A detailed analysis of the uniform system without a trap shows that the disordered states 
participate in a Bose glass phase. 
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Recent advances in experiments with ultracold atoms 
in magneto-optical traps have opened a new frontier in 
the study of strongly correlated systems. Some of the 
more remarkable early experimental breakthroughs in- 
clude the realization of Bose-Einstein condensation^] 
and fermionic superfluidity. More recent achievements 
include tuning a degenerate Fermi gas across BCS- 
BEC crossover via Feshbach resonance, and detect- 
ing superfluid-to-Mott insulator transitions of trapped 
bosons in optical lattices. 0, Q An unprecedented con- 
trol over experimental parameters in these systems makes 
them ideally suited for studying many-body phenomena. 

Until recently, the majority of experiments with 
trapped ultracold atoms have been performed on clean 
systems. Indeed, the ability to create perfect (defect free) 
optical lattices is a major advantage over condensed mat- 
ter experiments. On the other hand, the complete control 
over the trapping potentials makes it possible to intro- 
duce disorder and tune the disorder strength in a con- 
trolled fashion. The exciting possibilities of being able 
to study novel disorder-related phenomena such as An- 
derson localization and to explore novel quantum glassy 
phases have made the investigation of disorder in these 
systems an area of emerging interest. 

Disorder can be generated in optical lattices by ex- 
posure to speckle lasers liLZL adding an incommensu- 
rate lattice- forming laserQ or other means 01 • It 
is thus possible to investigate different disorder-induced 
phenomena in a controllable manner, in contrast e.g. to 
previous studies with granular superconductors or 4 He 
in vycor glass. The interplay between disorder and in- 
teractions in trapped Bose-Einstein condensates has re- 
cently been explored experimentally in 87 Rb, both in the 
continuum^, El 0] and in an optical lattice^]. On 
the theoretical front, such systems have been investi- 
gated within the frameworks of mean- field theories, 0,0, 
UUl,E| Bose-Fermi mapping. [lfij and the transfer matrix 
formalism. |l6j 

In this paper, we study interacting trapped bosons in 
one-dimensional optical lattices with a tunable random 




FIG. 1: Illustration of atoms in a potential trap. In the ab- 
sence of interactions between the atoms, disorder leads to an 
immediate localization of the single-particle wave functions. 



potential using a quantum Monte Carlo method. Previ- 
ous investigations have argued that under realistic con- 
ditions, such diagonal disorder dominates over random- 
ness in the hopping amplitude. || In experiments with 
optical lattices, the primary source of information about 
the state of the system arises from the analysis of the 
momentum distribution function, obtained from matter- 
wave interference after the release of the trap and subse- 
quent free evolution of the particles. In view of this, here 
we focus on the signatures of disorder on the momen- 
tum distribution. The simulations presented here pre- 
dict a unique enhancement of phase coherence at small 
to moderate strengths of disorder, when the groundstate 
has a Mott insulating region, i.e. a Mott plateau, at 
the center of the trap. Since Mott insulating states are 
now experimentally observable, and disorder strengths 
are controllable by tuning the intensity of the speckle 
laser, these results can be straightforwardly tested using 
currently available experimental methods. One techni- 
cal difficulty in using current setups is the typically long 
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wavelength of speckle patterns (typically 8-10 times the 
optical lattice spacing) and the limited size of optical lat- 
tices (40-65 lattice spacings along each axis). However, 
these difficulties are likely to be overcome in the near fu- 
ture with advances in realizing larger lattices and / or 
using more than one speckle laser to create disordered 
potentials with a shorter correlation. 

Bosons in an optical lattice are well described by the 
one-band Bose Hubbard model. 0] In the presence of a 
strong periodic (lattice) potential, single-particle Wan- 
nier functions localized on the lattice sites form a com- 
plete basis set. Interactions are typically not strong 
enough to excite higher vibrational bands, justifying the 
single band approximation. Hence the low-energy physics 
of trapped bosons in a one-dimensional optical lattice can 
be described by the Hamiltonian 
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Here b\(bi) creates (annihilates) a boson at site i, rij = 
b\bi is the number operator, and U is the strength of the 
repulsive on-site Hubbard interaction between bosons. 
The bare chemical potential controls the filling of the 
lattice, Vt is the strength of the trapping potential, a 
denotes the lattice spacing, and Wi introduces diagonal 
disorder in the form of a random site energy. The hop- 
ping amplitude t is set equal to unity in the simulations, 
and all other parameters are implicitly expressed in units 
of t. In the experiments, Wi is commonly realized by a 
speckle laser. Furthermore, recent measurements indi- 
cate significantly correlated disorder, with a correlation 
length that is typically several times longer than the op- 
tical lattice spacing, due to the diffraction-limited imag- 
ing of the speckles co ming from a diffusion plate onto 
the trapped atoms. fill Il2j In view of this, we consider 
finitely correlated disorder with a random potential Wi 
that is distributed as 



W l = 0, WiW v =Mi, v . (2) 

The overbar denotes averaging over disorder realizations, 
and A is disorder strength. The indices I are measured 
in units of 5a. 

In the absence of disorder, the Bose Hubbard model 
JU in a trap has been extensively studied. [lf| Due to the 
spatially varying trapping potential, the system is never 
in a uniform phase. Depending on the strength of the 
trapping potential, Vt, the on-site interaction, U, and 
the density of particles, n, the groundstate can belong 
to one of two classes. At sufficiently small U/t and total 
density, n, it consists of a supcrfluid (SF) domain ex- 
tending across the entire system. The site-occupation, 
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FIG. 2: Effects of disorder on the momentum distribution 
n(q) and on the real-space density profile n(x) in the weak 
(upper panels) and strong (lower panels) coupling limit. The 
density profiles are vertically offset to enhance clarity of pre- 
sentation. 



rii varies continuously from zero at the edges to a finite 
maximum value at the center of the trap. The momen- 
tum distribution, n(q), features a sharp peak at q = re- 
flecting strong phase coherence across the entire system. 
The full width at half maximum of the zero-momentum 
peak is inversely proportional to the correlation length £. 
In a thermodynamic SF £ diverges, whereas in the one- 
dimensional confined SF £ oc N and n(q — 0) oc ^N. 

On the other hand, at larger U jt and sufficiently large 
n, the groundstate contains one or more Mott insulating 
(MI) regions with integer rii, along with domains of SF 
at the edges of the trap and in between the MI regions (if 
there are more than one such regions). In the present 
study, we consider the simplest case of one MI plateau 
at the center of the trap and two SF domains at the 
edges - the local density profile consists of an extended 
region at the center with m = 1 and two regions with 
< m < 1 near the trap edges. In the presence of MI 
region(s), the range of coherence is greatly reduced, and 
the momentum distribution shows a weak peak at q = 0, 
and a subsequent slow decrease in n(q) as a function of 
q. Depending on whether there exists an MI region at 
the center of the trap, we find the effects of disorder can 
be markedly different. 

We use the stochastic series expansion quantum Monte 
Carlo method to simulate the disordered Bosc-Hubbard 
model on finite-sized lattices. 01 The density profiles, 
n(x), and the momentum disribution function, n(q), are 
measured within the grand canonical ensemble, using suf- 
ficiently large values of the inverse temperature, (3, in or- 
der to obtain groundstate properties. Disorder averages 
over 400-2000 realizations are performed, depending on 
the parameters. The lattice sizes are chosen to be suffi- 
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ccntly large to ensure that the local density vanishes at 
the edge of the trap. We found that a lattice size L = 120 
is necessary for large values of U (> 8t) at strong disorder 
(W/U > 1.5), whereas L = 80 was sufficient for the other 
parameter sets. Results are displayed for 40 bosons in a 
parabolic trap with Vra 2 = 0.015i for a range of values 
of U ft, ranging from the weak to strong coupling limit. 

Fig. [21 shows the consequences of a disordered site- 
potential on the momentum distribution function n(q) 
and the real-space density profile n(x) for two values of 
the on-site interaction, U ft. The upper panels show the 
results for U/t = 2, when the groundstate in the clean 
limit (W = 0) consists of a single SF domain extending 
acroos the trap. Adding disorder produces localization 
centers around which the wave functions get localized. 
This in turn reduces coherence and hence results in a 
smaller zero-momentum peak. At sufficiently weak dis- 
order strengths, these effects are found to be relatively 
small in finite systems, and the groundstate retains a 
predominantly SF character - the deviation of the mo- 
mentum distribution from its clean limit is exponentially 
small. However, with increasing disorder strength, the lo- 
calization effects strengthen, and hence the momentum 
distribution starts to deviate significantly, i.e. the peak 
at n(q — > 0) becomes shorter and broadened. 

On the other hand, for U/t = 8 (lower panels), the 
groundstate of the system, in the clean limit, has a Mott 
plateau at the center and the effects of disorder are con- 
siderably more interesting. The wave functions in the 
MI region in the absence of disorder are localized due 
to the strong Hubbard interactions. The onset of disor- 
der weakens the correlation between particles, and thus 
leads to a partial delocalization of the strongly localized 
Wannier orbitals at small to moderate strengths of dis- 
order. Phase coherence hence increases as the MI re- 
gion "melts", reflected in a stronger q = peak in n(q) 
with increasing disorder. For larger disorder potentials, 
the Anderson localization effects dominate, and the zero- 
momentum peak starts to decrease. 

FigEI shows the variation of the strength of the zero- 
momentum peak as a function of disorder for four differ- 
ent values of U/t, ranging from the weak to the strong 
coupling limit. For intermediate values of the interaction 
strength, the peak at n(q = 0) increases at small to inter- 
mediate disorder (reflecting the melting of the central MI 
region), reaches a maximum, and finally decreases con- 
tinuously for larger W. This is to be contrasted with the 
weak coupling limit, where n(q = 0) remains fairly un- 
changed for small W , and then decreases rapidly at larger 
disorder strengths. The exact position of the maximum 
and its value, n max {q = 0) depends on U/t, Vr/t, and 
n. The peak position shifts to larger disorder strengths 
with increasing Hubbard intearction strength, indicating 
a competition between Mott localization and Anderson 
localization that depends on the ratio of W vs. U. The 
maximum height, n max (q = 0) and its ratio to the clean 
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FIG. 3: Dependence of the zero-momentum peak height, 
n(q — 0), on the disorder strength W for different values 
of U /t. At U/t = 5, the groundstate consists of a single SF 
domain, whereas for the other values, there exist a finite MI 
region at the center of the trap. 



system value, no(q = 0), decrease rapidly with U/t. At 
U/t = 20, i.e. close to the hard-core limit, no visible 
increase with W is noticablc in n(q = 0). 

What is the nature of the state that results from the 
onset of disorder? Does the partial delocalization of the 
MI states result in a SF or a glassy phase? To an- 
swer these questions, we analyze the effects of disorder 
in the absence of a trapping potential. At Vp =0, the 
groundstate is in a pure phase, and can be conveniently 
characterized by measuring global observables, such as 
the compressibility and the stiffness. For shallow trap- 
ping potentials, the system is adiabatically connected to 
the periodic case, such that the domains of the different 
phases retain their global characteristics to a large ex- 
tent. This allows us to study the effects of disorder in 
the pure phases and extend the results to correspond- 
ing domains in the trapped system. Let us focus on the 
global compressibility, k — (3((n 2 ) — (n) 2 , and the stiff- 
ness, p s , to differentiate between various possible phases 
- SF, MI or Bose glass (BG). The stiffness gives the 
response of the system to a uniform twist in the spa- 
tial boundary as p s = d 2 E/d<f> 2 , under which the ki- 
netic energy term in the Hamiltonian is replaced by 
— t J2( e ~ l ^b\ +1 bi + h.c). In practice, the stiffness is sim- 
ply related to the fluctuations in the winding number in 
the simulations as p s = (W 2 ) /2/?L.[2£j The results shown 
in Fig. 0] are for uniform SF (left panels) and MI (right 
panels) phases in the presence of potential disorder. The 
dependence of n(q = 0) on W is qualitatively similar 
compared with the case of trapped atoms. For the SF 
phase, n(q = 0) remains practically unchanged at small 
W, and starts to decrease monotonically for larger disor- 
der. On the other hand, in the MI phase n(q = 0) first 
increases and then decreases with increasing W, confirm- 
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cause of this competition there is an intermediate regime 
with enhanced phase coherence. An analysis of the com- 
pressibility and stiffness in the corresponding periodic 
systems, i.e. in the absence of a trapping potential, re- 
veals that the resulting disordered state in either case 
is a Bose glass. The predicted enhancement of phase 
coherence by diagonal disorder suggests an interesting 
extension of recent experiments on optical traps in the 
presence of a tunable speckle laser. |l2| 
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FIG. 4: Dependence of the momentum distribution function 
at q = 0, the stiffness and the compressibility of SF (left 
panels) and MI (right panels) phases on disorder strength. For 
the SF phase, the groundstate is predominantly SF at small 
to moderate disorder and a Bose glass at large disorder. For 
the MI phase, the groundstate is a Bose glass at all strengths 
of disorder studied. The strength of the zero-momentum peak 
varies non-monotonically, but the groundstate always remains 
a Bose glass. 

ing that the disordered states are qualitatively similar to 
their counterparts in the presence of a trapping poten- 
tial. The stiffness and compressibility data allow us to 
further characterize the disordered states. For the SF 
groundstate, both p s and k are finite for the range of W 
over which n(q = 0) remains close to its W = value, 
consistent with a predominantly SF character of the dis- 
ordered state. At larger W, as n(q = 0) starts to deviate 
significantly, the compressibility remains finite, but the 
stiffness rapidly decreases to zero. Hence the groundstate 
at large disorder is a compressible insulator, i.e. a Bose 
glass (BG). From the present data it is not clear if one 
needs a finite critical disorder to destroy superfluidity or 
simply the localization length at small disorder is larger 
than the system sizes considered here. For the MI phase, 
both p s and k vanish in the clean system. With the onset 
of disorder, the stiffness remains zero, but the compress- 
ibility acquires a finite value. The disordered MI is thus 
a BG at all strengths of disordered considered here. 

In summary, based on quantum Monte Carlo simula- 
tions ensembles of interacting trapped atoms are found 
to respond in a highly non-trivial fashion to diagonal dis- 
order. If the clean system consists of a single superfluid 
domain, the particles are simply localized by the random 
potential. This is reflected by progressive suppression 
of the q — ► peak in the momentum distribution func- 
tion. On the other hand, if the clean system contains 
Mott insulating regions the response to on-site disorder 
is non-monotonic, and depends on the relative strength of 
the Hubbard interaction and the disorder potential. Be- 
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